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Abstract. 

We have conducted multi-epoch EVN observations of the 22.2 GHz water masers towards the high-mass young 
stellar object in AFGL 5142. With four observing epochs, spanning a time of ~1 year, 12 distinct maser features 
have been detected, 7 out of these detected in more than one epoch. The positions and the velocities of the 
VLBI features agree well with those of the emission centers previously identified by means of VLA observations. 
For a few features, persistent over three or four epochs, accurate values of the proper motions are derived. The 
observed proper motions have an amplitude of 15-20 km s~^, significantly larger than the range of variation of 
the line-of-sight velocities (±4 km s~^ around the systemic velocity). On the basis of their spatial distribution, 
the observed maser features can be divided into two groups. A model fit to the positions and velocities of the 
maser features of Group I, detected in the same region (within ~500 mas) where the massive YSO should be 
located, demonstrates that these might arise on the surface of a nearly edge-on Keplerian disk, rotating around a 
massive young stellar object. The maser features of Group II, found at large distances from the YSO (> 1"), have 
positions and line-of-sight velocities in agreement with the blue-shifted lobe of a large scale molecular outfiow 
(traced by the 1100"*" and SiO emission), and might result from the interaction between the gas fiowing away 
from the young stellar object and the ambient gas of the progenitor molecular core. 

Key words, masers - stars: formation - ISM: individual objects (AFGL 5142) - ISM: kinematics and dynamics 
- Radio lines: ISM 



1. Introduction species, such as OH, H2O, CH3OH, observed in the prox- 
5— i ■ imity of the high-mass proto-star. Multi-epoch VLBI ob- 
Ci ; The star formation process is better understood for low- gg^^ations can provide accurate (relative) positions, line- 
mass stars (~1 Mq) than for high-mass stars (> 10 Mq). transversal velocities of the maser spots (the 
The massive stars are less numerous, on average more individual, mas-scale compact, centers of emission), so 
distant from the Sun, and enter the ZAMS phase still ^^^^ ^j^^ 3-dimensional velocity distribution of the gas 
enshrouded m their progenitor dust and gas envelope, ^^^^^^ ^j^^ ^^^^^ transition can be derived. Moreover, 
making optical and near-mfrared observations imprac- yLgi measurements of the emission properties of the 
ticable. This explains why, to date, only a handful of ^-^^^^ ^^^^^ ^^^^^ ^^^^ ^g^j^Uy compared with the 
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maser excitation models to constrain the relevant phys- 
ical and geometrical parameters. 



the formation process of massive stars, high resolution The first multi-epoch VLBI experiments were per- 

observations at radio, millimetre and far-infrared wave- formed about 20 years ago towards a few of the 

lenghts are needed. The highest resolutions (< 1 mas) are strongest 22.2 GHz H2O interstellar masers in the 

obtainable through the Very Long Baseline Interferometry Galaxy, i.e. ^ ^ Orion-KL jCenzel et alj |l981a|) W5 1 

(VLBI) technique at radio wavelenghts, which can be dCenzel et al. 1981bl), Sagittarius B2 jReid et al.||l98a), 

used to observe the maser transitions of several molecular W49 (Gwinn et al..,1992,), and allowed both the determi- 

nation of the characteristic pattern of the kinematics of 

Send offprint requests to: C. Goddi, these regions and, by comparing the line-of-sight veloci- 

e-mail: cgoddi@ca.astro.it ties with the proper motions of the spots, the derivation 
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of accurate source distances. Since then, VLBI observa- 
tions of 22.2 GHz water masers have been mostly focused 
on those sources which, selected on the basis of interfero- 
metrical observations of (thermal) molecular tracers, are 
considered to be among the best high-mass protostellar 
candidates. The Very Long Baseline Array (V LBA) obser- 
vations carried out bv lMoscadelli et all 1)2000^) towards the 
water masers in IRAS 20126-1-4104, one of the best stud- 
ied examples of a high-mass Young Stellar Object (YSO) 
associated with a Keplerian disk and a jet/outflow system 
ijCesaroni et al.lll997) . suggest that the H2O masers could 
arise on a bipolar conical surface, excited by the interac- 
tion of an ionized jet with the surrounding molecular gas. 
VLBA observations of the 22.2 GHz masers performed 
towards intermediate-mass (NGG 2071 ISeth et al.ll2002() 
and high-m ass (IC 1396N. ISlvsh et alJll999l: W3 mS 5, 
llmai et all Eo 00 'I YSOs indicate that this maser emission 
can originate both in the protostellar disks and at the base 
of the molecular outflows. 

This paper presents multi-epoch VLBI observations 
of the 22.2 GHz H2O masers towards the high-mass star 
forming region A FGL 5142 (I RAS05274-H3345). At a dis- 
tance of 1.8 kpc JSnell et all FigSS). its far-infrared lumi- 
nosity is estimated to be 3.8 x 10'^ Lq ijCaroenter et alJ 
2002l(. Using the V ery Large Array (VLA) at 8.4 GHz 
Tofani et al T ill 99,^ bave revealed a faint (1 mJy), almost 
unresolved, continuum source (best interpreted as free- 
free emission from an ionized wind), which subsequently 
has been found to be coincident in position with the cen- 
ter of a CO bipolar outflow and with the origin of a iet 
obser ved in the H2 near-infrared emission (.Hunter et al.l 
19951). Owens Valley Radio Observatory (OVRO) data by 
Hunter et al.l l)l999^ show: 1) a well-coUimated SiO-jet 



and an HCO^-outflow, both aligned with the axis of the 
CO-outflow, and emanating from the 8.4 GHz continuum 
source; 2) a compact 88 GHz continuum source coincident 
in position (within the observational errors) with the 8.4 
GHz continuum. The radio flux and the bolometric lumi- 
nosity (estimated using the IRAS fluxes) can be explained 
if the exciting source is a massive object, with spectral 
typ e B2 or earlier. The VLA NHs-observations performed 
by IZhang et all 1(20021) show a 1" diameter (1800 AU) 
compact structure, whose morphology and kinematics are 
compatible with a rotating disk surrounding a high-mass 
young star. 

The 22.2 GHz water masers in AFGL 5142 have 
been observed with the VL A at three epochs (1991, 
TorrcUe s et all Il9 92!: 1992, iHunter et all Il995t 1998, 



Hunter et al 1 119991) ■ and found to be distributed within 
a few arcseconds from the 8.4 GHz continuum source. 
However, the VLA angular resolution ('^0."1) is inade- 
quate to determine the detailed spatial distribution of the 
maser spots, to measure their proper motions, and to in- 
vestigate the kinematics traced by the H2O masers in this 
source. At the distance of AFGL 5142, the linear resolu- 
tion attainable using the VLBI technique at 22.2 GHz is 
~2 AU, sufficient in principle to resolve and determine the 
velocity profile of a protostellar disk (of size of hundreds 




-5 

Velocity (km/s) 

Fig. 1. The Effelsberg total-power spectrum observed in 
Oct 1996 (connected plot) is compared with the integrated 
flux densities (histogram plot) of the VLBI channel maps. 
The velocity resolution of both spectra is equal to the 
channel width of 0.12 km s^^. 



of AU) and/or to study the mechanism of coUimation and 
acceleration of a jet. 

Section 2 of this paper describes our multi-epoch VLBI 
observations and gives technical details of the data anal- 
ysis, while Section 3 presents the observational results. 
Section 4 describes plausible kinematical models for inter- 
preting the measured positions and velocities of the maser 
features. The conclusions are drawn in Section 5. 

2. Observations and data reduction 

AFGL 5142 was observed in the 616 — 523 H2O maser 
line (rest frequency 22235.080 MHz) using the European 
VLBI Network (EVN) at four epochs (October 1996, and 
June, September, November 1997), each epoch consisting 
of 13 scans (of 6.5 minutes) distributed over an hour an- 
gle of 11 hours. The antennae involved in the observations 
were Medicina, Cambridge, Onsala, Effelsberg, Metsahovi, 
Noto, Jodrell and Shanghai^, with a subset of 5-7 of these 
(including always Effelsberg) observing at each epoch. 
Before and after each scan on AFGL 5142, several con- 
tinuum sources were observed for calibration purposes 
(0528-fl34, 0642-f449, 1803+784, 2251-^158, 3C454.3, 
OJ287). Both circular polarizations were recorded with 
a 1 MHz bandwidth centered at the LSR velocity of - 
4.8 km s~^. The data were processed with the MKIII cor- 
relator at the Max-Planck-Institut fiir Radioastronomie 
(Bonn, Germany), obtaining 112 spectral channels with a 
separation of 0.12 km s~^. 



^ The baselines to Shanghai have been almost completely 
flagged during the process of self-calibration. 
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Feature 


Vlsr 




Aa 


A<5 


14 


Vy 


T^mod 






(km s-i) 


(Jy) 


(mas) 


(mas) 


(km s-i) 


(km s-i) 


((km s-i) 




1 


-7.2 


16.9 


189.7 (0.3) 


1020.7 (0.2) 


23.7 (12.6) 1" 


-33.8 (8.2) 


41.2 (9.9) 




2 


-4.8 


5.9 


157.7 (0.2) 


1010.43 (0.03) 


4.2 (1.8) 


-15.2 (0.4) 


15.8 (0.6) 




3 


-4.8 


0.5 


157.3 (0.2) 


1012.1 (0.2) 








Group I 


4 


-3.6 


0.8 


159.0 (0.2) 


1017.26 (0.03) 










5 


-2.0 


4.1 


-32.5 (0.2) 


1290.15 (0.09) 










6 


-1.2 


1.0 


174.3 (0.3) 


1007.5 (0.2) 


4.4 (2.7) 


-19.1 (1.3) 


19.6 (1.4) 




7 


-0.6 


0.3 


28.7 (0.2) 


1251.03 (0.07) 










8 


0.7 


0.3 


55.5 (0.2) 


1245.85 (0.08) 










9 


-5.4 


0.3 


-50.3 (0.3) 


-18.99 (0.02) 








Group II 


10 


-5.2 


4.8 


0.00 


0.00 


0.00 


0.00 


0.00 




11 


-3.9 


3.0 


-655.2 (0.3) 


-566.7 (0.1) 


-25.3 (15.0) 


-9.8 (4.8) 


27.1 (14. 1) 




12 


-3.7 


1.2 


-179.3 (0.2) 


-75.6 (0.1) 


-1.0 (1.9) 


0.6 (1.3) 


1.2 (1.8) 



t The itaUcs indicate tentative values of proper motion components for features observed at only two epochs 



Data reduction was performed using the NRAO AIPS 
package, following the standard procedure for VLBI line 
data. Total power spectra of continuum calibrators were 
used to derive the bandpass response of each antenna. The 
amplitude calibration was performed using the informa- 
tion on the system temperature and the gain curve of each 
antenna. We did not use the "template spectrum" method 
(which consists of comparing total power spectra of differ- 
ent scans) owing to insufficient signal-to-noise ratio (SNR) 
of the total power spectra of the maser source. 

For each observing epoch, a single scan of a strong 
calibrator was used to derive the instrumental (time- 
independent) single-band delay and the phase offset be- 
tween the two polarizations. After removing the instru- 
mental errors, all calibrator scans were fringe-fitted to de- 
termine the residual (time-dependent) delay and the fringe 
rate. The corrections derived from calibrators were applied 
to the strong maser feature used as a phase reference. 

In the data analysis of each observing epoch we used 
as phase-reference the same maser feature at Vi^sR — ^ 
4.8 km s~^, which exhibits a simple spatial structure con- 
sisting of a single, almost unresolved spot. The visibilities 
of the phase-reference channel were fringe-fitted to find 
the residual fringe rate produced both by differences in 
atmospheric fluctuations between the calibrators and the 
source, and by errors in the model used at the correlator. 
After correcting for the residual fringe rate, the visibilities 
of the reference channel were self-calibrated to remove any 
possible effects of spatial structure. Finally, the corrections 
derived from the reference channel were applied to data 
of all spectral channels. 

Before producing synthesis aperture maps, we searched 
for maser emission in a wide field (12" x 12") centered 
on the reference feature using the fringe-rate mapping 
method (.WalkeniigSlI) . The detected emission centers, all 
within a few arcseconds from the phase center, were conve- 
niently mapped by producing (naturally weighted) maps 
extended over a sky area of {Aa cosd x Ad) 4" x 4" and cov- 



ering the velocity range from —10.5 to 0.7 km s^^. The 
CLEAN beam was an elliptical gaussian with a typical 
FWHM size of 2.1 x 1.1 mas. In each observing epoch, the 
RMS noise level on the channel maps, a, is close to the the- 
oretical thermal value, 0.03 Jy beam~^, for channels where 
no signal is detected, and increases to 0.3 Jy beam^^ for 
channels with the strongest components. 

Every channel map has been searched for emission 
above a conservative detection threshold (in the range 
5-10 cr), and the detected maser spots have been fit- 
ted with two-dimensional elliptical Gaussians, determin- 
ing position, flux density, and FWHM size of the emission. 
Hereafter, we use the term of "feature" to indicate a col- 
lection of spectrally and spatially contiguous maser spots. 
A maser feature is considered real if it is detected in at 
least three contiguous channels, with a position shift of the 
intensity peak from channel to channel smaller than the 
FWHM size. Fig. ^ compares the Effelsberg total-power 
spectrum observed in Oct 1996 with the integrated flux 
densities of the VLBI channel maps for the same epoch. 
One notes that, excepted for the velocity range from —7 
to —8 km s~^, the difference between the total-power flux 
and the flux recovered in the imaged VLBI field of view is 
within the amplitude calibration errors (± 50%). 

The relative positional uncertainty of the single maser 
spot is estimated using the expression 

AO = — FWHM, (1) 

where FWHM is the un-deconvolved spot size, / is the 
peak intensity an d a is the no-signal rms of the map 
l(R,eid et alJll98 8^. For most of the spots the positional 
uncertainty is of the order of ~50 - 100 /xas. 

3. Observational Results 

Counting all four epochs, 26 H2O maser features were de- 
tected. Several of these show a good agreement in rela- 
tive positions (within few mas) and line of sight velocities 
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Fig. 2. Measured relative proper motions of H2O maser features in AFGL 5142. For each of the time persistent features, 
the top and the bottom panels report the time variation of respectively the declination and the right ascension offsets 
(relative to the feature " 10"). In each panel, the dotted line shows the proper motion calculated by the (error- weighted) 
linear least-squares fit of the positional offsets with time. Among the features observed at three or more epochs, the 
proper motions are derived only for those moving in a straight line at constant velocity (within the positional errors). 
Adopting such a criterion, no proper motion is derived for feature "5". Tentative values of the proper motions are 
calculated also for the two features (label number "1" and "11") observed at only two epochs. 
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(within 1 km s^^) for two or more epochs, and therefore 
we assume that in these cases we identified features which 
persist over time. A final set of 12 distinct features, 7 out 
of these observed for more than one epoch, has been iden- 
tified. 

Table 1 gives the parameters of the features, deter- 
mined by fitting an elliptical Gaussian to the intensity dis- 
tribution of the maser spots contributing to the features' 
emission at different velocity channels. Col. 1 gives the fea- 
ture label number. Cols. 2 and 3 report respectively the 
line-of-sight velocity, Vlsr, and the integrated flux den- 
sity, Fint, of the highest-intensity channel, both averaged 
over the observational epochs for the time-persistent fea- 
tures. For each feature. Cols. 4 and 5 report the positional 
offset (of the first epoch of detection) calculated with re- 
spect to the feature with label number "10" , present at all 
four epochs. The positional offsets of a given feature are 
estimated from the (error-weighted) mean positions of the 
contributing maser spots. The bracketed numbers are the 
relative positional uncertainties, evaluated by taking the 
weighted standard deviation of the spot positions. With 
the 1-2 mas angular resolution of the EVN at 22 GHz, 
the emission of most of the spots is found to be spatially 
unresolved. Using a distance value of 1.8 kpc for AFGL 
5142, an upper limit to the spot size of ~2 AU is derived. 

The proper motions have been calculated performing 
a (error-weighted) linear least-squares fit of the positional 
offsets with time. Fig. |5| shows the time variation of the 
right ascension and declination offsets (relative to the fea- 
ture "10") for the persistent features. Among the features 
observed at three or more epochs, the proper motions are 
derived only for those moving in a straight line at con- 
stant velocity (identified with the label numbers "2", "6" 
and "12"). Tentative values of the proper motions are cal- 
culated also for the two features (with label number "1" 
and "11") observed at only two epochs. Cols. 6, 7 and 8 
of Table 1 report respectively the projected components 
along the R.A. and DEC axis, and the absolute value of 
the derived proper motions. The numbers in italics refer 
to features observed at only two epochs. The bracketed 
numbers are the formal errors of the linear least-squares 
fit. 

Fig. El compares our VLBI results with previous inter- 
ferometric observations. Top panels show the high- velocity 
molecular outflows seen in HC O"*" (1 — > 0) and SiO {v — 0, 
2^1) with OVRO ( Hunter et al.lll999li . The area com- 
prising almost the totality of the 22.2 GHz water maser 
emission detected by our VLBI observations, is indicated 
by a small filled rectangle at the center of the field of 
view. This area is expanded in the lower panel of Fig. 
which shows the spatial distribution of the VLBI maser 
features superimposed on top of the distribution of the 
VLA emission centers found bv lHimter et al.l 1 199,4 [l99^ 
at two different epochs. The measured proper motions are 
indicated by the arrows. 

The absolute position of the VLBI map has been de- 
termined by aligning the VLBI and VLA maser emission 
centers. The VLBI map has been shifted on top of the VLA 



map to find a shift that minimizes the root mean square 
difference between the positions of the VLBI and VLA 
spots. We find that each of the detected VLBI features 
has a good positional (within the VLA beam of ~100 mas) 
and line-of-sight velocity correspondence with one of the 
observed VLA spots. Conversely, we have detected VLBI 
features corresponding to each VLA emission center with 
peak flux density > 0.4 Jy. The good overlap found be- 
tween the VLBI and the VLA emission makes us confident 
that the absolute position derived for the VLBI map is ac- 
curate within the VLA positional uncertainty. 

On the basis of their spatial distribution, the VLBI 
maser features can be divided into two groups. Group I, 
comprising the first eight features, is found in the north- 
east corner of the area plotted in the lower panel of Fig. |3 
in the same region (within ~500 mas) where the 8.4 GHz 
and 88 GHz continuum emissions are detected. Group II 
includes the last four maser features more detached (> 
1") from the continuum emissions, extending towards the 
south and the southwest corner of the plotted area. 

Fig. 3 and Table 1 of iHunter et all l|l999lh show that 
at epoch 1998 two VLA emission centers (their compo- 
nents "6" and "7") are detected at positions separated 
about 3-4 from the 8.4 GHz continuum peak. We pro- 
duced VLBI maps also at these two locations, toward both 
directions detecting a single, weak (<Jy beam^^) spot in 
the last two epochs (September and November 97). No 
features persisted over time and, consequently, no proper 
motions are measured. The components " 6" and " 7" were 
among the s trongest ones in th e epoch 1998 VLA ob- 
servations of iHunter et al.l ljl999f) b ut were not detected 
in th e two prior VLA ep ochs H 991 . iTorrelles et aLl Frngi 
1992. lHunter et al.i l995V Our VLBI observations, extend- 
ing over the years 19 96 and 199 7, fall in between th e 1992 
I Hunter et aLlll99,'Tl) and 1998 l|Hunter et al.lll999l) VLA 
runs, and witness the first appearance of water maser 
emission in these two regions. Being so detached from the 
8.4 GHz and 88 GHz sources, the maser emission in these 
two locations very likely traces a site of star formation 
other than the one responsible for the line and continuum 
emission shown in Fig. 13 A single maser spot does not al- 
low us to derive information on the gas kinematics and we 
will not consider further these two regions in the following 
discussion. 



4. Discussion 

Fig. 13 shows that the spatial distribution of all the VLBI 
and VLA maser features is extended along a north— south 
direction. The outflow detected in t he CO, HCO"*" and 
SiO emission has a similar orientation (^H unter et al.ll995l 
■1999 ). In addition, although the line-of-sight velocities of 
the VLBI maser features do not vary smoothly, the mean 
velocity of the Group I cluster (toward the north) is more 
positive (redshifted) than that of the Group II (toward 
the south). Over an area of sky sligthly larger than that 
plotted in Fig. 13 a similar variation of the line-of-sight 
velocities (more redshifted (blueshifted) toward the North 
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Fig. 3. The upper panels show the contour maps of the high-velocity molecular line emission of HCO+ (1 0) {left) 
and SiO (w = 0, 2 ^ 1) (right), observed with OVRO in 1998 by Hunter et al. (1999; their Fig. 10). The area at the 
center of the OVRO field of view, where most of the 22.2 GHz water emission detected by our VLBI observations is 
concentrated, is expanded in the lower panel. This shows the s patial distribu tion of the VLBI features superimposed 
on top of the distribution of the VLA emission centers found bv lHunter et al.l l(l99f)Lll99i^ at the two different epochs 
(1992, open diamonds; 1998, filled dots). Each VLBI feature is identified with the label number given in Col. 1 of Table 
1. Different colours are used to distinguish the line of sight velocities of the features, according to the colour- velocity 
conversion code shown on the right-hand side of the panel. The arrows indicate the measured proper motions, with 
dashed lines used in case of more uncertain values. The dotted triangles drawn around the proper motion vectors 
represent the amplitude and orientation uncertainty. The region in the northeast corner enclosed by the dotted lines 
is shown in greater detail in the box at the center of the lower panel, whose coordinates are arcsec offsets relative to 
the position: RA = 05''27™30'; DEC = -h33°45'40". For the feature labeled "12", with a small measured value of the 
proper motion, cannot be assigned a definitive direction of motion. The amplitude scale for the proper motions is given 
at the bottom left corner of the panel. The small and the big cross es give respectively the position uncertainties of the 
VLA 8.4 GHz and the OVRO 88 GHz continuum emission peaks I Hunter et al.lll999|l . For each of the two Groups of 
EVN maser features, the black dot-dashed line indicates the axis whose average distance from the feature positions is 
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(South)) is also noteable in the VLA 22 GHz observa- 
tions of Hunter et al. (1999; their Table 1), where the 
most northward VLA maser compo nents (labeled " 9" , 
"10" and "11" in iHunter et~al]ll999ip occur at velocities 
(> -1 km s^^) higher than th e most southwarth compo- 
nents (labeled "2" and "8" in lHunter etaOliggal ) (< - 
3 km s~^). The observed variation of water maser velocity 
is qualitatively in agreement with the velocity distribution 
seen at a much larger angular scale in the HCO"*" and SiO 
maps (Fig. 131). A simple interpretation might be that all 
the detected (VLA and VLBI) maser features are tracing 
the flow motion in the innermost portion of the molecu- 
lar outflow. To make the discussion a bit more quantita- 
tive, considering that the line-of-sight velocity dispersion 
of the large-scale (diameter ~ 50") molecular outflow is 
'^lOO km s~^, and assuming a Hubble flow (velocity in- 
creasing linearly along the outflow axis) , one would derive 
a rate of line-of-sight velocity dispersion caused by the 
outflow of ^2 km s~^ arcsec"^. If that might explain the 
velocity dispersion of the Group II features (1.7 km s~^ 
over a distance of 1"; see Table 1), the Group I features 
show a much higher velocity dispersion (8 km s~^) across 
a smaller distance (0."35). 

The maser features of Group I have a sky-projected 
distance < 500-1000 AU from the 8.4 GHz and 88 GHz 
continuum sources, and should emerge near to the ex- 
pected location of the massive YSO. At such close distance 
from the YSO, it might be possible that the flow motion 
has not yet reached a stable configuration and turbulence 
might play a role in increasing the velocity dispersion of 
the gas. Alternatively, the maser features of Group I might 
move under the influence of the gravitational field of the 
massive YSO. Within a region of radius of ~1000 AU 
around the forming high-mass star, the current theory of 
star formation predicts that an accretion disk should be 
found. Looking at Fig. 13 one notes that the maser fea- 
tures of Group I have an elongated spatial distribution 
(the dot-dashed line indicates the elongation axis) and 
that the measured proper motions have orientation close 
to that of the elongation axis. This geometrical condition 
is what one would in principle expect if these maser fea- 
tures traced a rotating disk seen edge-on. However, one 
should note that the derived proper motions are relative 
to the feature "10" (not belonging to the Group I maser 
cluster), and, in order to obtain the absolute transversal 
velocities of the gas, one has to correct for the (unkwown) 
absolute motion of this feature. 

Even if only a small number of maser features is de- 
tected toward the Group I cluster, neverthless our ac- 
curate knowledge of their positions and line-of-sight ve- 
locities offers the chance to fit their motion using a 3-D 
Keplerian disk model. The model's free parameters are: 
the sky-projected coordinates of the YSO (at the disk cen- 
ter) ; the position angle and the inclination angle with the 
line-of-sight of the disk axis; the YSO mass. For a given 
set of input parameters, one can compute the position 
and the velocity vector of each feature and compare the 
model velocities with the observed velocities. The best fit 



to the data was obtained minimizing the ^ given by the 
squared sum of the error- weighted differences of the model 
and observed line-of-sight velocities. Incidentally we note 
that the same fit solution is found when the is cal- 
culated including also the two measured proper motions 
(for the features labelled "2" and "6" in Table 1). 

Looking for the disk axis orientation over the full 47r 
solid angle, the best fit solution is found with the disk 
seen almost edge-on (inclined 12° from the line-of-sight) 
and oriented on the sky parallel with the elongation axis 
of the Group I features (at P. A. = 153°). The fitted value 
of the YSO mass, My so = 38 Mq ± 20 M©, although 
determined with high uncertainty, strongly indicates that 
the central object is a massive YSO (M > 10 Mq). This 
result is in ag reement with that of previous observations 
llHunter et ab l ll 99,5lll 999h . indicating an exciting object 
of spectral type B2 or earlier, for whic h the theoretically 
expected value of the ma ss is > 10 Mq l|Vacca et al 
IPalla fc Zinneckeill2002l) . Adopting a distance of 1.8 kpc 
to AFGL 5142, the range of disk radii traced in our model 
by the maser emission extends from '-^30 AU to ^800 AU, 
which is consistent with the size of several hundreds of 
AU expected for an accreting disk around a massive YSO. 
From the barely resolved 88 GHz emission, interpreted in 
terms of optically th in thermal emission from a dusty core, 
iHunter et alJ l(l999l) derive a core size of ^5000 AU and 
a core mass of ^^145 Mq. These values are compatible 
with the results of our model, which tells us that within a 
radius of ^1000 AU the mass in Keplerian motion is < 

60 Mq. 

Recently IZhang et al ] 1I2OO2I) have observed the high- 
mass star- forming region AFGL 5142 in several NH3 inver- 
sion transitions using the VLA array. In corr espondence 
of the 88 GHz source of IHunter et all l)l999l) . they find 
a compact (1" in diameter), hot (70 K), NH3 structure, 
with a broad line emission, inter preted as an unres olved 
rotating disk. Looking at Fig. 2 of lZhang et a one 
sees that the velocity range (from -8 km s~^ to 2 km s~^) 
over which the NH3 "disk" emission is detected matches 
well with the velocity range of the Group I maser emis- 
sion. The value of th e "disk" mas s estim ated from the 
NH3 measurements of IZhang et al ] f2002') is 4 Mq, as- 
suming an NH3 abundance (relative to H2) of 10~^. This 
value of the "disk" mass, being much lower than the fitted 
YSO mass (38 Mg), is compatible with the Keplerian disk 
model proposed to explain the kinematics of the 22 GHz 
maser features that we observed. 

Owing to their large distance (2000-3000 AU) from the 
YSO believed to be responsible for the continuum emis- 
sions and the acceleration of the large scale molecular on- 
flow, the Group II of maser features might be also asso- 
ciated with a distinct (as yet undetected) YSO. However, 
the fact that the positions and the radial velocities of 
these features are in agreement with the blue-shifted lobe 
of the molecular outflow, makes us favour the intepreta- 
tion that their emission is excited by the interaction of 
the gas outflowing from the YSO with the ambient gas 
of the progenitor molecular core. Looking at Fig. 13 one 
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notes that the spatial distribution of the maser features of 
Group II is elongated (the dot-dashed line indicates the 
elongation axis) along a direction (at P. A. = 44°) that dif- 
fers from that of the large scale HCO+ and SiO outflows, 
oriented approximately north-south. Towards the infrared 
sour ces IRSl an d IRS3 in the NGC2071 star-forming re- 
gion, [Set^^^L] f|2002,l note similar differences of orienta- 
tion between the outflow structures on large (> 1000 AU) 
and small (10-100 AU) scales. This effect might be ex- 
plained thinking of either the density gradients in the am- 
bient medium, causing large-angle bends of the protostel- 
lar jets, or multiple, small scale outflows, whose merging 
creates the large scale flow. 

5. Conclusions 

Using the EVN we have observed the 22.2 GHz H2O 
masers towards the massive star forming region AFGL 
5142 for four epochs (from October 1996 to November 
1997). Previous high-angular resolution observations of 
several thermal tracers in the cm and mm-wavelength 
band indicate that the exciting source of the water maser 
emission is a high-mass YSO, of spectral type B2 or ear- 
lier. 

We identified a final set of 12 distinct water maser fea- 
tures, 7 out of these detected in more than one epoch. For 
a few features, persistent over three or four epochs, ac- 
curate values of the proper motions are derived. The ob- 
served proper motions have amplitudes of 15-20 km s^^, 
significantly larger than the range of variation of the line- 
of-sight velocities (±4 km s~^ around the systemic veloc- 
ity) . The positions and the velocities of the VLBI features 
agree well with those of the emission centers previously 
identified by means of VLA observations. 

On the basis of their spatial distribution, we have di- 
vided the maser features into two groups. Group I, com- 
prising eight VLBI features, is found in the same region 
(within ^500 mas) where the massive YSO, believed to 
be the exciting source of two compact, continuum sources 
(at 8.4 GHz and 88 GHz), should be located. A model fit 
to the positions and velocities of these features indicates 
that they might arise on the surface of a Keplerian disk 
seen nearly edge-on. The fitted value of the YSO mass, 
My so — 38 Mq ±20 Mq, although determined with high 
uncertainty, is in agreement with the results of previous 
observations. Maser features of Group II, found at larger 
distances (> 1") southward from the YSO, have positions 
and line-of-sight velocities in agreement with the blue- 
shifted lobe of a large scale molecular outflow (detected 
in the HCO+ and SiO emission), and might be excited by 
the interaction of the gas outflowing from the YSO with 
the ambient gas of the progenitor molecular core. 

This work demonstrates that multi-epoch EVN obser- 
vations are able to measure the proper motions of the 
strongest and longer-living 22.2 GHz maser features. Our 
EVN observations towards the source AFGL 5142 suf- 
fered two major drawbacks: 1) the time separation be- 
tween consecutive epochs (> 3 months) is too long com- 



pared to the average life time of the maser features; 2) 
modest sensitivity (with an average detection threshold of 
~0.35 Jy beam~^), owing to the fact that only 5-7 anten- 
nae (out of the 11 available to observe at 22.2 GHz) took 
part in each run. Next we plan to use the VLBA to better 
constrain the kinematical scenario suggested by the EVN 
observations, using a shorter time separation (~1 month) 
between two consecutive epochs and higher sensitivity. 
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